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Abstract 

We study by numerical simulations on a lattice the behaviour of the gauge-invariant two- 
point correlation functions of the gauge-field strengths across the deconfinement phase 
transition, both for the pure-gauge SU(3) theory and for full QCD with two flavours. 
Quenched data agree within errors with previous determinations, but have much higher 
statistics. A best-fit analysis of the data has been performed, both for the quenched 
and the full-QCD case, showing that the electric gluon condensate drops to zero at the 
deconfining phase transition. 
(PACS codes: 12.38. Gc, ll.10.Wx) 
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1. Introduction 

Gauge-invariant correlation functions of the field strenghts in the QCD vacuum play an 
important role in high-energy phenomenology and in stochastic models of QCD, both at 
zero temperature [1, 2, 3] and non-zero temperature [4, 5, 6]. For a recent review see Ref. 
[7]. Some years ago, a determination of such correlators at finite temperature was done 
on a 16 3 x 4 lattice, for the pure-gauge SU (3) theory, in a range of distances from 0.4 to 
1 fm approximately [8, 9] . The technique used to make the computation feasible is a local 
cooling of the configurations: this procedure freezes local fluctuations, leaving long-range 
correlations unchanged. In this paper, prompted by the progresses of the stochastic- 
vacuum approach to QCD [10], we improve the determination of the correlators at finite 
temperature T ~ T c for pure gauge SU(3), by use of a larger lattice (32 3 x 6) and bigger 
statistics. We also compute the correlators at T = and at T ~ T c in full QCD with 2 
staggered dynamical quarks. In Sect. 2 we recall the notation. The numerical results are 
presented in Sect. 3, while Sect. 4 contains an analysis of the data and a discussion. 



2. Notation 

To simulate the system at finite temperature, a lattice is used of spatial extent N a 3> 
N T , N T being the temporal extent, with periodic boundary conditions for gluons, and 
antiperiodic boundary conditions for fermions in the temporal direction. The temperature 
T corresponding to a given value of (5 = 6/g 2 is given by 



where a is the lattice spacing. In the quenched case a only depends on the coupling (3 
and, from renormalization group arguments, 




(2.1) 



(2.2) 
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where f{j3) is the so-called scaling function and is the scale parameter of QCD in the 
lattice regularization scheme. At large enough f3, f((3) is given by the usual two-loop 
expression: 

/8 \ 51/121 , 4 v 

fW) = exp (-— tt 2 /?) [1 + 0(1/ (3)] , (2.3) 

for gauge group SU(3) and in the absence of quarks. The expression (2.3) can also be 
used in a small enough interval of /3's lower than the asymptotic scaling region, and then 
Al is an effective scale depending on the position of the interval considered. For the range 
of values of /3's that we have considered (see below), its value, extracted from the string 
tension [11, 12], is about 4.9 MeV. 

The gauge-invariant two-point correlators of the field strengths in the QCD vacuum 
are defined as [1, 2, 3] 

T> mva {x) = g 2 (Tr[G, p (0)S(0,x)G ua (x)S\0,x)]) , (2.4) 



where G w = T a G^ p is the field-strength tensor, T a are the generators of the algebra of 
the gauge group in the fundamental representation, and 

S(0, x) = P exp (ig £ dt x^A^xt^j (2.5) 

is the Schwinger parallel transport from to a; along a straight-line path;* A p = T a A a ^. At 
zero temperature, that is on a symmetric lattice N a = N T , the correlators are expressed in 
terms of two independent invariant functions of x 2 , known as V(x 2 ) and V^x 2 ) [1, 2, 3]: 

T^n P ,ua{x) = (5 p J pa - 5 pa 5 pu ) [V(x 2 ) + Vx(x 2 ) 

dV 1 (x 2 ) 

~\~{X p Xyb ptj XpXfjbpy -\- XpXfjbpy XpXybp^a) ' (^'^) 

At finite temperature (N a 3> N T ) the 0(4) space-time symmetry is broken down to the 
spatial 0(3) symmetry and the bilocal correlators are now expressed in terms of five 
independent functions [4, 5, 6]. Two of them are needed to describe the electric-electric 

*Recently a strong dependence of the correlators on the shape of the path in the Schwinger string has 
been found numerically [13]. However, the stochastic-vacuum models naturally select the straight-line 
path. For that reason, we limit our present analysis to the straight-line Schwinger string. 
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correlations: 



g 2 {Tr[E i (x)S(x,y)E k (y)S\x,y)}) 

,dD? 



ik 



DE + D * + <-Sul 



+ UiU k 



dDf 

du 1 



(2.7) 



where E { = G, A is the electric field operator and u^ — x^ — [u 2 — (x — y) 2 ]. 
Two further functions are needed for the magnetic-magnetic correlations: 



g 2 (Tr[B i (x)S(x,y)B k (y)S\x,y)}) 



= 5, 



th 



D» + D? + u 



du 2 



UiU k - 



du 2 



(2.8) 



where B k = \eij k Gij is the magnetic field operator. 

Finally, one more function is necessary to describe the mixed electric-magnetic corre 
lations: 

1 dD* E 



g 2 (Tr[E i (x)S(x,y)B k (y)S\x,y)}) 



" n £ikn'U"n~ <^ 

2 0U4 



2.9) 



In Eqs. (2.7), (2.8) and (2.9), the five quantities D E , D E , D B , Df and D BE are all 
functions of u 2 , due to rotational invariance, and of u\, due to time-reversal invariance. 

From the conclusions of Refs. [4, 5, 6], one expects that D E is related to the (temporal) 
string tension and should have a drop just above the deconfinement critical temperature 
T c . In other words, D E is expected to be a kind of order parameter for confinement; 
on the contrary, Df does not contribute to the area law of the temporal Wilson loop. 
Similarly, D B is related to the spatial string tension [4, 5], while Df does not contribute 
to the area law of the spatial Wilson loop. 

The above arguments hold both for the quenched and the unquenched theory, with a 
suitable modification of Eq. (2.2) and (2.3). In particular for the full-QCD case, in order 
to fix the scale, we have used the lattice spacing as determined in [14]. 
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3. Results 



We have determined the following four quantities [15] 



V E (u 2 ,0) 



V E (u 2 ,0)+V E (u 2 ,0)+u 2 
V E (u 2 ,0)+V E (u 2 ,0) ; 



(u 2 ,0); 



du 2 



(3.1) 



Vf(u 2 ,0) 
V B (u 2 ,0) 



V B (u 2 ,0) +V B (u 2 ,0) +u 
V B (u 2 ,0)+V B (u 2 ,0) , 




(u 2 ,0) ; 



du 2 



(3.2) 



by measuring appropriate linear superpositions of the correlators (2.7) and (2.8) at equal 
times («4 = 0). Concerning the mixed electric-magnetic correlator of Eq. (2.9), it 
vanishes both at zero temperature and at finite temperature, when computed at equal 
times (w 4 = 0), as a consequence of the invariance of the theory under time reversal. 

We have chosen a 32 3 x 6 lattice (so that, in our notation, N T = 6) for the quenched 
case. The critical temperature T c for such a lattice corresponds to (3 C ~ 5.89 [16]. For full 
QCD we have used 2 species of staggered fermions with bare mass a ■ m = 0.0125 and 
a 32 3 x 8 lattice, for which T c corresponds to a coupling f3 c ~ 5.54 [17]. The standard 
"R-version" of the HMC algorithm [18] has been used in the full-QCD case. 

For the quenched theory the behaviour of T> E and T> E is shown in Figs. 1 and 2 
respectively, at different values of T/T c with the physical distance in the range from 
~ 0.25 fm up to ~ 1.25 fm. A clear drop is observed for V E and T>f across the phase 
transition, as expected. The analogous behaviour for V B and Vf is shown in Figs. 3 
and 4. No dramatic change is visible across the transition for the magnetic correlations. 
These behaviours of the correlators were already known from Refs. [8, 9], with which we 
agree within the errors. In Figs. 1-4 the thick continuum line has been obtained using 
the parameters of the best fit at T = obtained in Ref. [8], Eqs. (2.10) and (2.11). 

In Figs. 5-8 we present the analogous data for the case with dynamical fermions. In 
full QCD previous results at T = have been reported only for the case of 4 flavours [19]. 
Therefore, in order to compare with T = 0, we have performed a simulation at f3 = 5.55 on 
a 16 4 lattice: the results are reported in the same figures. Also in this case it is apparent 
that T> E and T> E stay almost constant at their zero-temperature value up to the phase 
transition, where they undergo a sharp drop. Instead no dramatic change is visible for 
the magnetic correlators. 
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Our results, both for the quenched and the full-QCD case, are in agreement with those 
already found in Refs. [8, 9] and can be summarized as follows: 

(1) In the confined phase (T < T c ), up to temperatures very near to T c , the correlators, 
both the electric-electric type (2.7) and the magnetic-magnetic type (2.8), are nearly 
equal to the correlators at zero temperature: in other words, D E ~ D B ~ D and 
Df ~ Df ~ D 1 for T < T c . 

(2) Immediately above T c , the electric-electric correlators (2.7) have a clear drop, while 
the magnetic-magnetic correlators (2.8) stay almost unchanged, or show a slight 
increase. 

In the next section we shall report on a quantitative analysis of the data displayed in 
Figs. 1-8. 



4. Quantitative analysis 

Inspired by our previous analyses of the correlators in the T = case (see Refs. [15, 8, 19] 
and also [20]), we have performed best fits to the lattice data for the correlators (3.1) and 
(3.2) at finite temperature T (and at equal times, i.e., u± — 0), with the functions (here 
x = \u\): 

V E (x)=A e-^ x + ^ , vf(x) = A 1 e-^ x + ^ , (4.1) 
a; 4 

V B [x) = Boe-^ + — A , Pf (x) = B ie -^ BX + — A , (4.2) 

X 4 X 4 

where, of course, all the coefficients must be considered as functions of the physical tem- 
perature T. The four independent functions (4.1) and (4.2) are written as the sum of a 
non-perturbative exponential term plus a perturbative-like term behaving as l/x A (in fact, 
a term of this form is predicted by ordinary perturbation theory).* In Refs. [8, 19] the 
perturbative-like term had the form ~ e~^ aX /x 4 ": the exponential term e~^ aX has been 

*Of course the coefficient of l/x 4 is regularization-scheme dependent. In Eqs. (4.1) and (4.2) we refer 
to the lattice regularization scheme; other schemes could give different values [21]. 
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neglected here (i.e., we fix /i a — 0), since, in the spirit of the Operator Product Expansion 
(OPE), we will concentrate on the behaviour of the correlators at short distances. 

Two quantities of physical interest enter in our best fits to the lattice data: 

(1) The correlation length of the gluon field strengths, defined as A^ = I/ha for the 
electric correlators and \b = 1//ib for the magnetic correlators. (A priori, we 
distinguish the two correlation lengths A^ and A# in the parametrization (4.1)- 
(4.2). However, as we shall be below, we have found that our data can be well 
fitted using the same correlation length A^ = A# for the electric and the magnetic 
correlators.) 

(2) The gluon condensate, defined as 

G 2 = : G a G a :) (a. = £) . (4.3) 



These two quantities play an important role in phenomenology. The correlation length is 
relevant for the description of vacuum models [1, 2, 3]. 

The gluon condensate was first introduced in Ref . [22] , in the context of the S VZ sum 
rules. As pointed out in Refs. [19, 20], lattice provides us with a regularized determination 
of the correlators. We shall briefly repeat here the argumentation originally reported in 
Refs. [19, 20]. As in Ref. [22], our correlators can be given an OPE [23]: 

2^W(*) x Zo Cl{x){1) + C ^ X) ° 2 + £ C A x ) m A : Wf ■■) + • • • > ( 4 - 4 ) 

if we have Nf quark flavours with masses rrif. (Of course, the last term in Eq. (4.4), i.e., 
the mixing to the quark condensates, is absent in the quenched theory.) The mixing to the 
identity operator C\{x) has al/i 4 behaviour at small x. The mixings to the operators of 
dimension four C g (x) and Cf(x) are expected to behave as constants for x — > 0. Higher- 
order terms in the OPE (4.4) are neglected. The coefficients of the Wilson expansion are 
usually determined in perturbation theory and are known to be plagued by the so-called 
"infrared renormalons" (see for example Ref. [24] and references therein). In the same 
spirit of Ref. [22], we shall disregard the renormalon ambiguity, as we did in Ref. [19]. 
With the normalization of Eq. (4.4), this gives C g (0) ~ 1. The contribution from the 
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quark operators in (4.4) can be neglected because it is higher order in 1/(3. The gluon 
condensate is then, using the parametrization (4.1)-(4.2) for the correlators: 

G 2 = ^(A + A 1 + B + B 1 ) . (4.5) 

G 2 is the sum of an electric contribution, that we shall call G 2 ele \ plus a magnetic contri- 
bution, that we shall call G 2 ma9 ^ , which at non-zero temperature T are in general different 
and should be distinguished: 

G (ele) _ Tr[ ^ 2] . G (rna 9 ) = tf, * ^ 

IT 2 TV 2 

When using the parametrization (4.1)-(4.2) for the correlators, one easily finds that: 

G 2 ele) = - 2 {A, + A{) ; Gt a9) = ~ 2 {B + B 1 ) . (4.7) 

Let us discuss now the results obtained from the best fits to our data with the functions 

(4.1) -(4.2). First of all, we have tried a best fit to the data for the magnetic correlators 

(3.2) with the functions (4.2), where the mass fis of the non-perturbative exponential 
terms has been put equal to the corresponding value obtained in the T = case: 

Hb = 4.53(7) fm _1 (quenched) ; 

Hb = 3.5(2) fin" 1 (full QCD) . (4.8) 

The quenched value has been taken from Ref. [8], Eqs. (2.10) and (2.11); the correlators 
at T = are reproduced by the thick lines in Figs. 1-4. Instead, the full-QCD value 
has been extracted from a best fit to the new data at T — (f3 — 5.55 on a 16 4 lattice), 
reported in Figs. 5-8. The magnetic correlators are well fitted, up to distances of about 
x ~ 0.8 fm, by the functions (4.2) with the mass [Ib fixed to the value given in Eq. (4.8). 

The results obtained for all the various cases are reported in Table I. From these 
results one sees that all the coefficients, B , Bi, b and b±, are rather stable, when vary- 
ing the temperature T: only the non-perturbative coefficient B Q shows a slight increase 
when increasing T, so that, by virtue of Eq. (4.7), we can say that the magnetic gluon 
condensate G^ ma9 * ) slightly increases across the transition at T c [see Fig. 9]. Let us also 
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observe that the pure-gauge magnetic gluon condensate is sensibly higher than the cor- 
responding value in the full-QCD case: this is in agreement with the fact that the gluon 
condensate G 2 is expected to increase with the quark mass [25], tending towards the 
asymptotic (pure-gauge) value, as already checked on data at T = [19]. The stability of 
the perturbative coefficients b and b\ is expected since they are UV-dominated cut-off 
terms. More precisely, these perturbative coefficients are practically independent on the 
lattice volume (N% x N T ), as one can verify by direct computation [26]. In general, (at a 
fixed lattice volume N% x N T ) they depend on f3 and it is known that bo = Oil/0 2 ) and 
b\ = 0(1/ /3). However, in the range of values of T that we have considered, A/3//3 <C 1 
and the /3 dependence of the perturbative coefficients can be neglected within the errors. 
For the same reasons, these magnetic perturbative coefficients should be practically equal 
to the corresponding coefficients a and a\ in the electric correlators (4.1). 

Therefore, on the basis of these considerations, we have tried a best fit to the data for 
the electric correlators (3.1) for distances from 3 up to 5-6 lattice spacings (corresponding 
approximately to the range of physical distances 0.3 -j- 0.6 fin) J with the functions (4.1), 
where the perturbative coefficients a and a x have been fixed to the (weighted) average 
values of the corresponding magnetic coefficients b and b\ reported in Table I for all the 
temperatures that we have considered: 

a = 0.55(2) ; a\ = 0.35(1) (quenched) ; 

a = 0.64(2) ; a x = 0.33(1) (full QCD) . (4.9) 

Moreover, we have fixed the mass ha of the non-perturbative exponential terms in (4.1) 
to the same value /xg used for the magnetic correlators [Eq. (4.8)], which is in turn the 
value at T = 0. In Table II we report the results obtained for the perpendicular electric 
correlator Vf [see Eq. (3.1)]. The coefficient A +Ai of the non-perturbative part of the 
correlator, which, by virtue of Eq. (4.7), is proportional to the electric gluon condensate 
G^ le \ sharply decreases across the transition [see Fig. 10]. Again, we find that the pure- 
gauge electric gluon condensate is sensibly higher than the corresponding value in the 
full-QCD case, in agreement with the general claim done in Ref. [25]. 

Finally, we have performed a best fit to the values for the difference 

Vf(x)-Vf(x) = --^L(x) (4.10) 
tThis is the range of distances where we have data for the parallel electric correlator at all temperatures. 
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between the two quantities displayed in Figs. 2 and 1 respectively (Figs. 6 and 5 for the 
full-QCD case). The results are reported in Table III. Since (from previous experience 
at T = 0) the quantity in the r.h.s. of Eq. (4.10) is expected to be dominated by the 
perturbative term, in this case we have left also the coefficient free parameter in 

the fit, in order to test the validity of the assumption that the perturbative coefficients 
are temperature independent. The test is perfect for the quenched case and reasonable 
for the full-QCD case. One finds that the coefficient Ai of the non-perturbative part of 
the function Vf stays compatible with zero, within the (large) errors, across the phase 
transition at T c . So the clear drop seen in the quantities T> E and T>f_ across T c seems 
to be entirely due to the coefficient A of the non-perturbative part of the function 
V E alone. This result is consistent with Refs. [8, 9] and again confirms the conclusion 
of Refs. [4, 5, 6], where V E (or, better, its non-perturbative part) was related to the 
(temporal) string tension <je- It was also shown in Refs. [4, 5] that V B (or, better, its 
non-perturbative part) is related to the spatial string tension a s . Existing lattice results 
[27, 28, 16] indicate that a s is almost constant around T c and increases for T > 2T C : this 
fact is in good agreement with the behaviour that we find for the coefficient B of the 
non-perturbative part of the function V B (see Table I). 

In summary, our best-fit analysis of the data leads to the following conclusions: 

(1) The correlation lengths = 1/ha for the electric correlators [Eq. (4.1)] and A# = 
1/fiB for the magnetic correlators [Eq. (4.2)] are equal and do not change across 
the deconfining phase transition at T c . 

(2) The electric gluon condensate, defined in Eqs. (4.6) and (4.7), drops to zero at 
T = T c , whilst the magnetic gluon condensate is practically unchanged, showing a 
small increase [see Figs. 9 and 10]. 

(3) As expected, the coefficients of the perturbative terms are temperature independent 
and are the same for the magnetic and the electric correlators. 
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TABLE CAPTIONS 



Tab. I. Results obtained from a best fit to the data of the magnetic correlators (3.2) with 
the functions (4.2), for the various temperatures that we have examined: "q" stands 
for "quenched" data, while "f" stands for "full-QCD" data. The mass ^ of the 
non-perturbative exponential terms has been put equal to the corresponding value 
obtained in the T = case [Eq. (4.8)]. An asterisk (*) near the value of some 
parameter means that the parameter was an input for the best-fit. 



Tab. II. Results obtained from a best fit to the data of the perpendicular electric correlator 
Vf [see Eq. (3.1)] with the functions (4.1), where the perturbative coefficients 
a and a\ have been fixed to the (weighted) average values of the corresponding 
magnetic coefficients bo and b\ reported in Table I for all the temperatures that we 
have considered [see Eq. (4.9)]. Moreover, we have fixed the mass fj, A of the non- 
perturbative exponential terms in (4.1) to the same value used for the magnetic 
correlators [Eq. (4.8)], i.e., to the corresponding value obtained in the T = case. 
The notation used is the same as in Table I. 



Tab. III. Results obtained from a best fit to the data of the difference between the perpen- 
dicular electric correlator Vf and the parallel electric correlator [see Eqs. (4.10) 
and (3.1)] with the functions (4.1), where the mass \ia of the non-perturbative ex- 
ponential terms in (4.1) has been fixed to the same value used in Table II. The 
notation used is the same as in Tables I and II. 
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Table I 



theory 


T/T c 


B 

(10 8 MeV 4 ) 


B l 

(10 8 MeV 4 ) 


(far 1 ) 


bo 


h 


X 2 /N 


q 


0.952 


1320(70) 


62(70) 


4.53 (*) 


0.54(5) 


0.28(3) 


0.7 


q 


0.974 


1290(73) 


-50(70) 


4.53 (*) 


0.60(5) 


0.36(3) 


1.2 


q 


1.007 


1240(52) 


-51(54) 


4.53 (*) 


0.57(4) 


0.36(3) 


0.6 


q 


1.030 


1436(63) 


24(77) 


4.53 (*) 


0.51(5) 


0.34(3) 


0.5 


q 


1.065 


1305(54) 


-55(70) 


4.53 (*) 


0.56(4) 


0.37(3) 


1.1 


q 


1.127 


1455(49) 


91(67) 


4.53 (*) 


0.50(4) 


0.35(3) 


1.9 


q 


1.261 


1490(38) 


17(60) 


4.53 (*) 


0.54(3) 


0.36(2) 


0.5 


f 


0.73 


446(25) 


37(16) 


3.5 (*) 


0.59(5) 


0.19(4) 


1.7 


f 


0.94 


461(23) 


24(24) 


3.5 (*) 


0.66(5) 


0.30(3) 


1.4 


f 


1.02 


463(18) 


-22(21) 


3.5 (*) 


0.63(4) 


0.37(2) 


0.63 


f 


1.18 


510(30) 


-11(27) 


3.5 (*) 


0.64(3) 


0.40(3) 


0.23 


f 


1.48 


574(35) 


8(34) 


3.5 (*) 


0.65(3) 


0.38(2) 


1.6 



Table II 



theory 


T/T c 


A + A t 
(10 8 MeV 4 ) 


Ha 
(frrT 1 ) 


a + cli 


X 2 /N 


q 


0.952 


1193(30) 


4.53 (*) 


0.90 (*) 


0.2 


q 


0.974 


930(25) 


4.53 (*) 


0.90 (*) 


1.6 


q 


1.007 


596(31) 


4.53 (*) 


0.90 (*) 


2.1 


q 


1.030 


318(42) 


4.53 (*) 


0.90 (*) 


0.9 


q 


1.065 


197(27) 


4.53 (*) 


0.90 (*) 


0.2 


q 


1.127 


56(60) 


4.53 (*) 


0.90 (*) 


0.3 


q 


1.261 


-88(57) 


4.53 (*) 


0.90 (*) 


0.4 


f 


0.73 


381(10) 


3.5 (*) 


0.97 (*) 


0.2 


f 


0.94 


413(20) 


3.5 (*) 


0.97 (*) 


0.84 


f 


1.02 


288(25) 


3.5 (*) 


0.97 (*) 


1.9 


f 


1.18 


186(30) 


3.5 (*) 


0.97 (*) 


3.4 


f 


1.48 


43(25) 


3.5 (*) 


0.97 (*) 


3.3 
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Table III 



theory 


T/T c 


A l 

(10 s MeV 4 ) 


(fnT 1 ) 




X 2 /N 


q 


0.952 


-70(305) 


4.53 (*) 


0.34(9) 


3 x 10" 3 


q 


0.974 


243(285) 


4.53 (*) 


0.34(8) 


1 x 10~ 3 


q 


1.007 


508(332) 


4.53 (*) 


0.36(7) 


2 x 10~ 2 


q 


1.030 


570(587) 


4.53 (*) 


0.37(15) 


1.4 x 10" 2 


q 


1.065 


544(790) 


4.53 (*) 


0.38(15) 


6 x 10~ 2 


q 


1.127 


502(829) 


4.53 (*) 


0.39(15) 


1 x lO" 1 


q 


1.261 


23(1733) 


4.53 (*) 


0.46(25) 


2 x 10~ 2 


f 


0.73 


29(61) 


3.5 (*) 


0.20(7) 


N = 


f 


0.94 


5(70) 


3.5 (*) 


0.33(7) 


0.01 


f 


1.02 


153(80) 


3.5 (*) 


0.37(7) 


0.1 


f 


1.18 


-10(40) 


3.5 (*) 


0.53(6) 


1.9 


f 


1.48 


-110(80) 


3.5 (*) 


0.55(6) 


1.2 
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FIGURE CAPTIONS 



Fig. 1. The quantity Vf [Eq. (3.1)], in units of MeV 4 , versus the physical distance (in fm), 
for different values of T/T c in the quenched theory. The thick continuum line has 
been obtained using the parameters of the best fit to the data at T = [Ref. [8], 
Eqs. (2.10) and (2.11)]. 

Fig. 2. The quantity Vf [Eq. (3.1)], in units of MeV 4 , versus the physical distance (in fm), 
for different values of T/T c in the quenched theory. The thick continuum line has 
been obtained using the parameters of the best fit to the data at T = [Ref. [8], 
Eqs. (2.10) and (2.11)]. 

Fig. 3. The quantity Vf [Eq. (3.2)], in units of MeV 4 , versus the physical distance (in fm), 
for different values of T/T c in the quenched theory. The thick continuum line has 
been obtained using the parameters of the best fit to the data at T = [Ref. [8], 
Eqs. (2.10) and (2.11)]. 

Fig. 4. The quantity Vf [Eq. (3.2)], in units of MeV 4 , versus the physical distance (in fm), 
for different values of T/T c in the quenched theory. The thick continuum line has 
been obtained using the parameters of the best fit to the data at T = [Ref. [8], 
Eqs. (2.10) and (2.11)]. 

Fig. 5. The quantity Vf [Eq. (3.1)], in units of MeV 4 , versus the physical distance (in fm), 
for different values of T/T c and for T = in full QCD. 

Fig. 6. The quantity Vf [Eq. (3.1)], in units of MeV 4 , versus the physical distance (in fm), 
for different values of T/T c and for T = in full QCD. 

Fig. 7. The quantity Vf [Eq. (3.2)], in units of MeV 4 , versus the physical distance (in fm), 
for different values of T/T c and for T = in full QCD. 
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Fig. 8. The quantity Vf [Eq. (3.2)], in units of MeV 4 , versus the physical distance (in fm), 
for different values of T/T c and for T = in full QCD. 

Fig. 9. The magnetic gluon condensate (^""^(T) [see Eqs. (4.6) and (4.7)], in units of 

versus T/T c . The black circles refer to the quenched case, while the 
white circles refer to the full-QCD case. 

Fig. 10. The electric gluon condensate G^iT) [see Eqs. (4.6) and (4.7)], in units of 
Q<fiie)(rp _ versus T/T c . The notation is the same as in Fig. 9. 
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